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Abstract 



We give reasons in support of the use of an effective intensity-dependent pion-nucleon 



in ' couphng Hamiltonian for describing the properties of the pion multiphcity distri- 

r-| ■ bution and the corresponding factorial moments within the thermal-density matrix 



approach. We explain the appearance of the negative-binomial (NB) distribution for 



Oh! 
(D . 

^ . pions and the well-known empirical relation of Wroblewski . Our model Hamiltonian 

• J • is written as a linear combination of the generators of the SU{1, 1) group. We find 

X 

. the generating function for the pion multiplicity distribution at finite temperature 

T and discuss the properties of the second-order factorial moment. Also, we show 
that an intensity-dependent pion-nucleon coupling generates the squeezed states of 
the pion field. At T = 0, these squeezed states become an inherent property of the 
NB distribution. 
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1 Introduction 



Recently a considerable amount of experimental information has been accumulated 
on multiplicity distributions of charged particles produced in pp and pp collisions 
in the center- of-mass energy range from lOGeV to ISOOGey . The Koba-Nielsen- 
Olesen (KNO) scaling [1], which was previously observed in the ISR cm. energy range 
from 11 to 63 GeV [2], was shown to be violated in the regime of several hundred 
GeV [3]. The violation of the KNO scaling is characterized by an enhancement of 
high-multiplicity events leading to a broadening of the multiplicity distribution with 
energy. 

The shape of the multiplicity distribution may be described either by its C 
moments, Cg — or by its central moments (higher-order dispersions), 

Dq^ ((n - g = 2, 3 . . .. The exact KNO scaling imphes that all Cq mo- 

ments are energy independent and that all central moments Dg satisfy a generalized 
Wroblewski relation [4,5] as (n) —>■ oo : 

Dq = Aq{n)-Bq, (1) 

with the energy-independent coefficients Aq and Bq. For the pp and pp inelastic 
data below lOOGeV, the coefficients Aq and Bq are approximately equal within 
errors. This fact implies that the elementary Poisson distribution resulting from the 
independent emission of particles is ruled out. 

The total multiplicity distribution Pn of charged particles for a wide range of 
energies (22 — 900GeV) is found to be well described by a negative-binomial (NB) 
distribution [3,6] that belongs to a large class of compound Poisson distributions [7]. 
The NB distribution is a two-step process [8] with two free parameters: the average 
number of charged particles (n) and the parameter k that affects the shape (width) 
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of the distribution. The parameter k is related to the dispersion D — D2hy the 
relation 



The observed broadening of the normalized multiplicity distribution with increasing 
energy implies a decrease of the parameter k with energy. The KNO scaling requires 
constant k. 

Although the NB distribution gives information on the structure of correla- 
tion functions in multiparticle production, the question still remains whether its 
clan-structure interpretation [8] is simply a new parametrization of the data or 
has a deeper physical insight [9]. Measurements of multiplicity distributions in pp 
collisions at TeV energies [10] have recently shown that their shape is clearly dif- 
ferent from that of the NB distribution. The distributions display the so-called 
medium-multiplicity "shoulder" [11,12]. A satisfactory explanation of this effect is 
still lacking [12]. 

In this paper we propose to study another approach to multiplicity distri- 
butions based on a pion-field thermal-density operator given in terms of an effec- 
tive intensity- dependent pion-nucleon coupling Hamiltonian with SU{\, 1) dynami- 
cal symmetry. We assume that the system of produced hadronic matter (pions) is 
in thermal equilibrium at the temperature T immediately after the collision. 

The paper is organized as follows. In Sec. 2 we present the basic ideas of our 
model. In Sec.3 we discuss the shape of multiplicity distribution, the correlations, 
and the Wroblewski relation. Finally, in Sec. 4 we conclude with a few remarks on 
the squeezing properties of the model in connection with the possible extension to 
include the two-pion coupling in the effective pion-nucleon Hamiltonian. 
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2 Thermal-density operator for the pion field 



At present accelerator energies the number of secondary particles (mostly pions) pro- 
duced in hadron-hadron collisions is large enough, so that the statistical approach 
to particle production becomes reasonable . Most of the properties of pions pro- 
duced in high-energy hadron-hadron collisions can be expressed simply in terms of 
a pion-field density operator. We neglect difficulties associated with isospin and, for 
simplicity, consider the production of single-mode pions. As a consequence of this 
restriction we are only able to calculate multiplicity distributions and multiplicity 
correlations of pions. The energy dependence of the multiplicity distribution will 
reveal the scaling properties of the collision dynamics, and the deviation of the dis- 
tribution from a Poisson distribution will reveal correlations between the produced 
pions. 

We expect that in high-energy collisions most of the pions are produced in the 
central region | y |< F, where Y = ln{s/w?) is the relative rapidity of the colliding 
particles. In this region the energy - momentum conservation has a minor effect if 
the transverse momenta of the pions are limited by the dynamics. 

The density operator poT for a free pion field in a heat bath of temperature T 

is 



1 



(3) 



Pot — 



where 



InZ — 



uj{a}a -\- A), 
-(5\uj-ln{l-e-^'^). 



(4) 



The quantity Aa; in Hq denotes the vacuum energy of the free pion system. The 



zero-point energy " corresponds to A = |. In the limit as T — > 0, the density 
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operator poT reduces to po =| 0)(0 | and represents the density operator for the 
pion-field vacuum state. 

The mean number of thermal (chaotic) pions is 

= (5) 

Owing to the interaction with the nucleon field the density operator poT is 
transformed by means of the unitary ^-matrix operator into 

Pt = SporS^ (6) 
1 

Z' 



where 

H = SHoSl (7) 

We regard H as an effective Hamiltonian describing the pion-nucleon system. At 
T — 0, the Pt becomes p — S \ 0)(0 | S^. The knowledge of pr gives the possi- 
bility of finding the pion-multiplicity distribution PT{n) at finite temperature T by 
calculating 

Prin) = (n I PT I n), (8) 
where | n) — (n!)~^/^a'''" | 0) are the pion number states. 

The first-order moment of Prin) gives the average multiplicity 

{n) = J2nPTin) (9) 

and higher-order moments of Prin) give information on dynamical fluctuations from 
(n) and also on multipion correlations. All higher-order moments can be obtained 
from the pion-generating function 

Gt(z) = Yl = Tr{pTZ^} (10) 
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by a suitable differentiation over where N — a^a. Thus the normahzed factorial 



moments are 



^ {n)i ^ ' dzi ^ ^ 



and the normalized cumulant moments K„ are 



K, - inr'^. (12) 
These moments are related to each other by the formula 

^. = E('^/)^.-^^^ (13) 

For the Poisson distribution, all the normalized factorial moments are identically 
equal to 1 and all cumulants vanish for g > 1. 

As a measure of dynamical fluctuations we shall mostly consider the 
q — 2 moments. They are directly related to the dispersion D of the multiplicity 
distribution Prin): 

F,-l^K, = {^f-j^ (14) 

The quantity F^ is also known as the two-particle correlation function g^*^^) = 
F2. Its numerical values are usually used to characterize the shape of the multiplicity 
distribution as well as to indicate the tendency of particles to bunch. The g^'^'' can 
take the following values: 



6 



gr(2) — 1 (^Poisson) 

{ 1 {sub — Poisson/antibunching) 

) 1 {super — Poisson/ bunching) (15) 

= 2 {chaotic) 

) 2 {enhanced bunching). 

The validity of the Wroblewski relation means that g^"^^ takes a value between 1 and 
2. 

The standard pion-nucleon interaction Hamiltonian is linear in the pion field 
operators. In our simplified model it corresponds to the Hamiltonian 

H ^uj{a''a + h+g{a'' + a), (16) 

where lo and g are real and positive parameters. This Hamiltonian can be written 
in the form (7) using the following S- matrix operator: 

^ = Dt(i) = exp{^(a-at)} (17) 

and A = | — (^)^ ■ The density operator pT corresponding to (16) describes displaced 
thermal states, i.e., the superposition of coherent and thermal 
states of pions. 

The generating function Gt{z) is [13] 

Gt{z) = — ^ — exp{-nc — }, (18) 

' l + (l-^)nT ^ '^l + {l-z)nT^' ^ ' 

where fic — {g/^^)^ denotes the average number of coherent pions. The normalized 

factorial moments Fq are given in terms of Laguerre polynomials: 

F, = g!(l + 7)-%(-7), (19) 
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where ^ — nc/riT- We note that the average number of pions is 

(n) ^UT + fic (20) 

and the two-pion correlation function is 

^(^) = 2 - {^f. (21) 

According to the relation (14), the square of the dispersion L>, in the limit (n) — > oo 
and 7 fixed, is 

(7^r-l-(T^r- (22) 
(n) 1 + 7 

It is always smaller than one as long as 7 7^ 0. 

On the other hand, by keeping the coupling constant g fixed and letting the 
temperature approach either zero or infinity ( variable 7 ), we find that g^'^'> — 1 
for coherent production and 51^^) = 2 for chaotic production. In this case, the 
Wroblewski relation (1) with a constant slope A cannot be obtained . To remedy 
this, we require that 7 is an energy-independent parameter [14]. This requirement 
means that the square of the coupling constant g'^ 

g'^u'-^in) (23) 
i + 7 

increases linearly with (n) or with T, which has the same effect. In this way, the 
Wroblewski relation can be satisfied for all energies. In analogy with quantum optics 
[15], we shall call the couphng of the type (23) intensity dependent. 

We propose to study the following form of the effective pion-nucleon intensity- 
dependent coupling Hamiltonian: 



H^e{N + X) + K{a\/N + 2A - 1 + h.c), (24) 

where e'^ = u"^ + Ak^ and A = (0 | ^ | 0)/e is related to the vacuum energy of the 
system. In the limit A — > 00 and k ^ 0, so that 2n^\ — g is, finite, the interaction 
part of the Hamiltonian H reduces to the standard pion-nucleon Hamiltonian (16). 



We observe that the operators 

Ko = N + X, (25) 



K_ = ayN + 2X-l, 



K+ = \/7V + 2A-la^ 



form the standard Holstein-Primakoff [16] reahzations of the su{l, 1) Lie algebra, 
the Casimir operator of which is 

C ^Kl-]-[K+K_ + K_K+\^X{X-l)i. (26) 

The Hamiltonian E. is thus a hnear combination of the generators of the 5'C/(1, 1) 
group: 

H^eK^^K{K^^K_). (27) 
The corresponding S-matrix that diagonahzes this Hamiltonian is 

^(^) =e-^(-^+-^-), (28) 

with 

t/i^ = ^. (29) 
The initial-state vector for the pion field, S{ff) | 0) =| 6*), is [17] 

I Q) = {l-th^e)>^U-m''{^^^^^^)'^' I n). (30) 

We note that the matrix element squared | {n \ 9) p can be written in the form of 
the NB distribution 



Pn^W = \{n\9)\' (31) 



r(n + 2A), n{e) 2A 



n!r(2A) 'n(^) + 2A' 'n(^) + 2A 
where n{e) ^ 2Xsh?{e). 



The full pion-density operator px = Pt(^) is 

PT{0)-T.PT{n)\n,e){n,ei (32) 

n 

where prin) — (riT)"/(l -\- utY^^ represents the Bose-Einstein (geometric) distri- 
bution function and \ n,9) = S{9) \ n). It is easy to see that the states | n, 9) form 
a complete orthonormal set of eigenvectors of the Hamiltonian H, i.e., 

H \n,e) ^uj{n + \)\n,e). (33) 

The calculation of the matrix element (n | S{9) \ m) — {n \ 9) is much 
faciliated if we present the operator S{9) in the antinormal form 

S{9) = e*'^W^- (c/i(^))2^°e-*'^W^+. (34) 

For n > m, we find 

(n I m,9) = {ch{9))-'\-th{9)r--{ ^^^^^^ (35) 
where P^'^^ are Jacobi polynomials [18]. A similar expression is obtained for n <m: 
{n I m,9) = {ch{9)r\th{9)r^^^!^^±^ (36) 
These matrix elements squared are used to find the pion multiplicity distribution 

PT{n)=J2PT{n)\{n\m,9)\^ (37) 

m 

and its generating function 

Gt(z) = 5^pT(m)(m,^ I 2* I m,^) (38) 

m 

— '^pT{m)z'^ I {n I 171,9) 1^ . 

m,n 
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3 Pion-generating function and its moments 



The average multiplicity (n), the dispersion and all higher-order moments (n'^), 
g = 1, 2, . . ., at the temperature T are obtained by a suitable differentiation over z 
from the pion-generating function Gt{z). The close analytic form of Gt{z) can be 
found by observing that 

(m, e\z^ \m,e)^{e\z^ \ e)y^Pj^'^^-'\x), (39) 

where 

^ z + (i-zrsh\e)ch'{e) 

z-(i-zysh^{e)ch^{e)' ^ ' 

z-{i- z)sh\e) 
^ " i + {i-z)sh^{e)' 

{e\z^\e) = Go{z) (41) 

m 

is exactly the generating function of the NB distribution with a constant shape 
parameter 2 A and the average number of pions equal to n{9). 

Taking into account the generating function of Jacobi polynomials [18] 

j2p(o,2X-i)^^-jym ^ 2''^-^R-\l + y + Ry-^\ (43) 



It is easy to see that 



R = sjl-2xy + y\ 

we obtain the final form of the pion- multiplicity generating function Gt{z): 

Gt{z) = G^{z){l + nT)-^2''^-^R^\l + + RtY'''^, (44) 
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where now yx — ^^T^ y. Using this generating function the following average number 
of pions and the multiplicity dispersion at the temperature T are found: 

(n) = n{0)+nT + \n{0)nT, (45) 

A 

A A 

where 

d\e) = ^n\d) + n{d). (46) 

Combining the two relations in (45), we can rewrite in the form of a 
quadratic function of (n): 

= A,{T) + A,{T){n) + A^{T){n)\ (47) 

with temperature-dependent coefficients ^j(T), i = 0, 1,2. In order to explain the 
Wroblewski relation , the coefficient A2{T), 

^(^- if)-' + if)- ■ (-) 

should be less than one. If A2{T) is fixed, say, by experiment, then it gives a 
definite relation between the vacuum energy A of the pion field and the number 
fiT of the produced thermal pions. Thus, for example,the data [9] on negatively 
charged particles produced in pp coUisions suggest A2{T) = 1/3. This value of the 
coefficient A2{T) restricts the average number of thermal pions to < 1 for all 
values of A > 3/2. 

We notice that at T = the Gt{z) becomes the generating function of the NB 
distribution, Gq{z), and 

L>2 (fie) 11 

+ (49) 



(n)2 T=o n^{e) 2A n{e) ' 
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as it is to be expected from the NB distribution . The shape parameter 2A of 
our NB distribution is given by the vacuum energy of the pion field in the nucleon 
environment, and has nothing to do with either the number of pion sources (cells) 
or the number of clans. The Wroblewski relation at T = 0, 

d(e) An(e) + B, n{e) > 1, (50) 

has energy-independent coefficients A = (2A)"^/^ and B = (A/2)^/^. If A > 1/2, we 
have A <1. 

For the temperature T going to infinity, we obtain 

= 2-(i + !f))- (51) 

r->-oo A 
= l + th^{26). 



in) 



This result shows that at very high temperature the distribution of pions will become 
chaotic if 9 is very small. 

The validity of the Wroblewski relation imphes that the produced pions have 
a tendency to bunch. In the standard approach, the bunching of pions is usually 
attributed to the presence of a quadratic (two-pion) interaction term in the Hamil- 
tonian (16), which is now of the form 

H = uj{a^a + ^)+ 9{a^ + a) + ^^(a^^ + a^). (52) 

Recently the properties of the density matrix pT, corresponding to this Hamiltonian, 

were investigated for the photon field [19]. It was found that px then defined the 
squeezed coherent thermal states [20]. 

Our effective intensity- dependent coupling Hamiltonian (24) is, however, 

highly nonlinear and nonquadratic in the pion field. It is , therefore, to be expected 

that it also generates squeezing that, in our model, should depend on the value of 
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the parameter A. Namely, we know that for A — oo the pion distribution reduces 
to the superposition of the coherent and thermal states that show no squeezing . 

We study the squeezing properties of the pion field in two Hermitian quadra- 
ture operators ai and a2 defined by 

a = ai + ia2, (53) 

which satisfy [01,02] = i/2. The corresponding uncertainty relation is AaiAa2 > 
1/4, where variances Aai^2 are defined by (Aai^2)^ = (c^i 2) ~ (^^1,2)^- A state of the 
pion field is considered squeezed if either Aoi or Aa2 are smaller than 1/2. If we 
define the relative variance with respect to (Aai 2)^0/1 = 1/4 as 

5i,2 = 4(Aai,2)' - 1, (54) 

then the squeezing condition becomes 

- 1 < < 0, i = l or 2. (55) 

At T = our model gives Q 

51 = 2((ata) + (at2) -2(at)2), 

52 = 2((ata) - (at2)), (56) 

where 

(at) = -th{e)J2p^^ie)V^^T2\, (57) 

n 

(at2) = th\e) + 2A)(n + 2A + 1). 

n 

^The finite-T case is considered elsewhere. 
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It is easy to see that for A » 1/2, the relative variance is Si > and. jS'2 ~ 0. 
However, for moderate values of A > 1/2, and such that 

^{n + 2A) {n + 2X + l)^(n + 2X) + l- + • • •, 

we find squeezing for 5*2 < 0. 
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4 Conclusions 



In this paper we have proposed to study an intensity- dependent pion-nucleon cou- 
pUng Hamiltonian with SU{1, 1) dynamical symmetry, within a multipion-produc- 
tion model in which the pion field is represented by the thermal-density operator. 

We have shown that this Hamiltonian explains in a natural way the appearance 
of the NB multiphcity distribution for pions and the Wroblewski relation .The shape 
parameter of the NB distribution is related to the vacuum energy of the pion field 
in the nucleon environment. 

Also, we have shown that, depending on the value of the parameter A, an 
intensity-dependent pion-nucleon coupling is able to generate the squeezed states 
of the pion field. At T = 0, these squeezed states | 9) then become an inherent 
property of the NB distribution. 

For T 7^ 0, we have found the explicit analytic form of the pion-multiplicity 
generating function that may be used for obtaining all higher-order moments of the 
pion field. 
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